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We study the time evolution of the entropy of a collapsing spherical domain wall, from the point 
of view of an asymptotic observer, by investigating the entropy of the entire system (i.e. domain 
wall and radiation) and induced radiation alone during the collapse. By taking the difference, we 
find the entropy of the collapsing domain wall, since this is the object which will form a black hole. 
We find that for large values of time (times larger than t/R s ~ 8), the entropy of the collapsing 
domain wall is a constant, which is of the same order as the Bekenstein-Hawking entropy. 
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I. INTRODUCTION 

Hawking's 1974 paper (see Refs.[l|, 0]) sparked great 
interest in both the existence of the radiation given off 
during the loss of mass for a pre-existing static black 
hole and speculations as to what his result implied for 
the validity of quantum field theory in the presence of a 
black hole. Hawking showed that the radiation given off 
was in fact thermal, which only increased the interest. 

In 1972 Bckcnstein argued that a black hole of mass 
M has an entropy proportional to its area (see Ref.[l3[). 
Further calculations by Gibbons and Hawking showed 
that the entropy of a black hole is always a constant, 
regardless of the type of metric which is used (see Ref. [|| ) . 
They showed that the expression for the entropy is given 
by 
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where A^or is the surface area of the horizon and R s is 
the Schwarzschild radius. 

Recently, subsequent work has been done in theories 
of quantum gravity to reproduce this result (see for ex- 
ample Refs.jil, [f|). So far, most of these calculations 
have involved only pre-existing black holes. The typi- 
cal method for finding the entropy of the black hole is to 
first determine the temperature of the black hole using 
the Bogolyubov method. Here, one considers that the 
system starts in an asymptotically flat metric, then the 
system evolves to a new asymptotically fiat metric. By 
matching the coefficients between the two asymptotically 
fiat spaces at the beginning and end of the gravitational 
collapse, the mismatch of these two vacua gives the num- 
ber of particles produced during the collapse. Hence, 
what happens in between is beyond the scope of the Bo- 
golyubov method. Therefore the time-evolution of the 
thermodynamic properties of the collapse cannot be in- 
vestigated in the context of the Bogolyubov method. 

In order to study the thermodynamic processes of the 
system, we will first develop the partition function. To 
do so we shall employ the so-called Liouville-von Neu- 
mann approach (see Refs. [Til. [l2|). This approach uti- 
lizes the invariant operator approach developed by Lewis 
and Reisenfeld (see Ref. jToj ) . In the invariant operator 
formalism, it was observed that any Hermitian operator 



I(t) satisfying the Liouville-von Neumann equation, 
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can be used to find the exact wavefunctions of the 
Schrodingcr equation, up to an overall time-dependent 
phase factor. In the Liouville-von Neumann approach, 
the time-dependent Hamiltonian is replaced by the time- 
dependent invariant operator, since the eigenvalues of 
this operator arc time-independent and all the time- 
dependence of the system is placed into a non-linear aux- 
iliary equation. 

To find the parameter /3, we shall use the so- 
called Functional Schrodingcr formalism (see for exam- 
ple Refs.@,H). In general, the Functional Schrodingcr 
formalism yields a functional differential equation for the 
wavefunctional, ^[g^v, X 11 , <f>, O], where g^ is the met- 
ric, X M the position of the object, $ is a scalar field and 
O denotes all the observer's degrees of freedom. Since 
the Functional Schrodingcr formalism depends on the ob- 
server's degrees of freedom, one can introduce the "ob- 
server" time into the quantum mechanical processes in 
the form of the Schrodingcr equation. The wavefunc- 
tional is then dependent on the chosen observer's time, 
hence one can view the quantum mechanical processes 
of a given system under any choice of space-time folia- 
tion. One benefit of this formalism is that one can solve 
the time-dependent wavefunctional equation exactly us- 
ing the invariant operator method, in particular in the 
present paper. Here we study the propagation of in- 
duced radiation, represented as a scalar field in the back- 
ground of gravitational collapse, which is governed by 
the harmonic oscillator equation with a time-dependent 
frequency. We solve the equations of motion to find the 
time-dependent wavefunctional for the system. One can 
then define the occupation number of the induced ra- 
diation using the wavefunctional, which is the gaussian 
overlap between the initial vacuum state and the wave- 
function at any given later time. By taking the occupa- 
tion number to be of the form of the Planck distribution, 
we can then define (3 which is also time-dependent. 

Since the partition function is time-dependent, one 
can then study the time-dependent thermodynamic pro- 
cesses associated with gravitational collapse. Therefore 
the Functional Schrodingcr formalism and Liouville-von 
Neumann approach together go beyond the approxima- 
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tions of the Bogolyubov method, since the system is al- 
lowed to evolve over time, which allows one to investigate 
the intermediate regime during the collapse. 

In this paper we will investigate the collapse of a spher- 
ically symmetric infinitely thin domain wall (represent- 
ing a shell of matter) . The details about the collapse will 
depend on the particular foliation of space-time used to 
study the system. For the purposes of the question we 
are studying, we are outside asymptotic observers. Thus, 
we study the collapse of the domain wall from the view 
point of a stationary asymptotic observer. 



II. MODEL 

To study a realization of the formation of a black hole, 
we consider a spherical Nambu-Goto domain wall (repre- 
senting a domain wall of matter) that is collapsing. This 
will be done using the Functional Schrodinger equation, 
where we will consider a minisuperspace version of the 
Wheeler-de Witt equation. Since we are interested in 
the entropy of the collapsing domain wall from the point 
of view of an asymptotic observer, the relevant degree of 
freedom for the collapsing spherical domain wall is the 
radial degree of freedom R(t). The metric of the system 
is then chosen to be the solution to Einstein's equations 
for a spherical domain wall. Therefore, outside the wall 
the metric is Schwarzschild, as follows from the spherical 
symmetry 

ds 2 =-fl-*k} dt 2 + fl - ^ 1 dr 2 + r 2 dQ 2 (3) 

where R s = 2GM is the Schwarzschild radius in terms of 
the mass M of the wall and 



By taking the ratio of these equations we then have 
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where R = dR/dt. 

By integrating the equation of motion for a spherically 
symmetric domain wall, Ipscr and Sikivic (sec Ref.p|) 
found that the mass is a constant of motion and is given 
by 



M = AiraR 2 



v/1 + R 2 - 2naGR 
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where a is the surface tension of the domain wall and 
R T = dR/dr. Before we proceed we wish to discuss the 
physical relevance of Eq. (TlT)|) . First consider the case 
where R T = 0, i.e. for a static domain wall. The first 
term in the square bracket is just the total rest mass 
of the domain wall. When the domain wall is moving, 
i.e. R T 0, the first term in the square bracket takes 
the kinetic energy of the domain wall into account. The 
last term in the square bracket is the self-gravity, or the 
binding energy of the domain wall. Therefore we can 
identify Eq. (jTUJ) as the total energy of the system, hence 
the Hamiltonian of the system. Thus, we will refer to 
Eq. (fTT)|) as the Hamiltonian. 

Using Eq.([7]) the Hamiltonian for the wall, Eq. lfTT))) , 
can be written as (see Ref.Q) 



H = 4tt<jB 3/2 R 2 
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where R = dR/dt. The canonical momentum near the 
horizon, the R ~ R s regime, is given by 



dn 2 = d9 2 + sin 
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By Birkhoff's theorem, the metric inside the domain wall 
must be flat, hence Minkowski 



ds z = -dT' z + dr 2 + r 2 dn 2 



(5) 



where T is the interior time. The interior time is related 
to the asymptotic observer time t via the proper time r 
of the domain wall. 
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and 
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where [i = er(l — 2ir<rGR). In this region the wall Hamil- 
tonian Eq. (|ll[) in terms of the canonical momentum is 
then 

infiB^R 2 

ia all ~ i 

VB 2 - R 2 



= ^(Bn R ) 2 +B(47rfiR 2 ) 2 (13) 

The Hamiltonian has the form of the energy of a rel- 
ativistic particle, \/p 2 + m 2 , with a position-dependent 
mass. 

In the region R ~ R s , Eq.jHJ) tells us that B ~ 0. From 
Eq.fTTS]) one can see that the mass term can be neglected, 
thus Eq. ([T3]) reduces to 



H W all ~ —BUr 



(14) 



where we have chosen the negative sign since the domain 
wall is collapsing. 
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Considering only the classical solution, the near- 
horizon solution of Eq. (fTTj) can be written as 



R w —B 
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Solving Eq. (|15p we then have 

R(t) ^R s + (R Q - R s )e- t ' R ' 



(15) 



(16) 



Eq. p^|) implies that a collapsing domain wall crosses its 
own Schwarzschild radius only after an infinite amount 
of asymptotic observer time t. 



III. RADIATION 

Here we consider the radiation given off during gravi- 
tational collapse. We will consider the radiation given off 
by the entire system and the particles produced during 
collapse alone. 

To investigate the radiation, we consider a scalar field 
<f> in the background of the collapsing domain wall. The 
scalar field is decomposed into a complete set of basis 
functions denoted by {/&(?")} 



$ = 5^o fc (t)A(r) 



(17) 



The exact form of the functions fk(r) will not be im- 
portant for us. We will, however, be interested in the 
wavefunction for the mode coefficients {a^}. 

The Hamiltonian for the scalar field modes is found by 
inserting Eq. (fTT)) and Eq.Q and Eq.© into the action 



-g\g^d^d^. 
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The Hamiltonian for the scalar field modes, arrived at 
from Eq. (fT8"|) . takes the form of uncoupled simple har- 
monic oscillators with independent mass and couplings 
due to the non-trivial metric. Using the principle axis 
transformation, the Hamiltonian can be diagonalized and 
written in terms of eigenmodes denoted by b. ft was 
shown in Ref. Q that in the regime of R ~ R s , the Hamil- 
tonian for a single mode can be written as 
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where is the momentum conjugate to b, and where m 
and K are constants whose precise values are not impor- 
tant to us. 



A. Entire System 



From Eq. p4|) and Eq. (jT9|) we can write the Hamilto- 
nian of the entire system as 



H = H wa u + Hb — —BUji + B- 
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where = —id/dR, Hb = —id/db. The wavefunction 
for the entire system is then a function of b, i?, and t, 
which we can write as 



* = R,t). 



(21) 



Using Eq. ([20"|) we can then write the Functional 
Schrodinger equation as 
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To solve Eq. (|2"2"|) we choose the classical background for 
the collapsing domain wall. Since the distance of the 
domain wall only depends on time, see Eq. (|15p . we can 
then write 
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dR~dt ~ 2^i~W + ~2 
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Subtracting the first time from both sides and grouping 
like terms, we can then write Eq. (|2"2"]) as 



b a 2 * k 
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Making use of the classical equation of motion, Eq. (|T5]) . 
Eq. ([23]l becomes 

b a 2 * 



2m db 2 



— b 2 ^ = 2i — . 

2 dt 



(24) 



Thus we treat the background classically and the radia- 
tion quantum mechanically. 

Dividing Eq. (|24p through by B, we now rewrite in the 
standard form of a harmonic oscillator 
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and 
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(26) 



(27) 



where Wg = K/m. Here we have chosen to set rj(t = 0) 
0. 

The solution to Eq. (|25[) is given by (see Ref.@) 



1/4 



exp 
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where p n = dp/dr] is the derivative of the function p(j]) 
with respect to 77, and p is the real solution to the non- 
linear auxiliary equation 



P A 

and where a is the phase, which is given by 

drf 



1 n 
«(»;) = -2 J 



(29) 



(30) 
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B. Radiation Only 

In this section we only consider the particles created 
during the collapse of the domain wall. To describe 
this we use only Eq. (fH))) . In this case the Functional 
Schrodingcr equation becomes 



where 



b a 2 * a' a* 

I h v[r = 7 

2m db 2 2 dt' 



(31) 



Dividing through by B and writing in the standard form 
yields 

i a 2 * 

dfj 



m , . ,<9* 



(32) 



where 



and 
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(33) 



(34) 



where ojq is defined as before. The solution to Eq. (j32"l) is 
the same as Eg . ([28)) except with the replacement r\ now 
given by 77 (Eq.([3"3"]l). i.e. 



1/4 



exp 



2 U p 2 



\np\ 

(35) 

Hence each wavefunction evolves with a different time 
parameter. 

C. Taking the Temperature 

For an observer, the complete information about the 
radiation, in the background of the collapsing domain 
wall, is contained in the wavefunction. Consider an ob- 
server with detectors that are designed to register parti- 
cles of different frequencies for the free field $. Such an 
observer will interpret the wave function of a given mode 
at some late time in terms of simple harmonic oscillator 
states. 

For brevity, we consider the following analysis using 
the notation for the entire system. However, the analysis 
for the radiation alone is the same. 

The wave function Eq. (f2"5)l (Eg. (|3"5|) ) can be decom- 
posed into suitably chosen vacuum basis wave functions, 
denoted by {4> n } at the final frequency Q. The number 
of quanta in eigenmode b can be evaluated by decom- 
posing the wavefunction in terms of the states {(p n } and 
by evaluating the occupation number of that mode. We 
decompose the wavefunction at a time t > tf in terms of 
a simple harmonic oscillator basis at t = 0. The decom- 
position is 



(36) 



(37) 



which is the overlap of a Gaussian with the t = sim- 
ple harmonic oscillator basis functions. The expectation 
value of the occupation number at given eigenfrequency 
u> is then given by 
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(38) 



Since the initial state is chosen to be the vacuum state, 
the number of quanta in eigenmode measured by the oc- 
cupation number is a consequence of only the collapse, 
since the observer will initially measure zero quanta. 

As stated previously, the observer will interpret the 
wavefunction of a given mode in terms of simple har- 
monic oscillator states. Thus the basis functions <f> n are 
chosen to be simple harmonic oscillator basis states at a 
frequency ui 



Mb) = 
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(39) 



where TL n are the Hermit c polynomials. Therefore using 
Eq. (|37p and Eq. (|39p we can write the inner product as in 
Ref. [7|, [8| , for n = odd c n = and for n = even 




(wp 2 )V4 



where 



P= 1- 



P 
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Substituting Eq.(|40j) and Eq.(|4Tj) and performing the 
sum over n = even in Eq. (|38p we then have 
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By fitting the occupation number to that of the Planck 
distribution, 



1 



(43) 



the occupation at eigenmode b can then be used to find 
the temperature of the radiation. From Eq. ([28|) and 
Eq. (f3"5]) we see that the occupation number for each of 
the systems will be different, since both ui and 77 are dif- 
ferent. 



IV. ENTROPY FUNCTION 

Here we develop the entropy function for studying the 
time evolution of the collapsing domain wall. Again for 
brevity we use the notation for the entire system, where 
again the analysis is the same for the radiation only sys- 
tem. 
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A. Partition Function 

To study the entropy of the system, we will first de- 
velop the partition function for the system. Following the 
procedure used in Ref . fl2| for the Liouville-von Neumann 
approach, we can write the partition function as 



Z = Tr [e~ 01 ] 



(44) 



where / is any operator which satisfies the Liouville-von 
Neumann equation Eq.Q and (3 is a free parameter. 
From Ref.0], we can write the invariant operator / as 



- + (nbP - mp v bf 



(45) 



We can therefore write the partition function as 
,1 



Z = Tr exp 



- + (^bP - mpr,b) 2 



(46) 



We note that we can rewrite the invariant as 

1 = U/i) Vp) -h™-™*® 
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where 



and 



n{t) = a){t)a(t) 
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Here n(t) is the time-dependent number of states. At a 
particular time t, one has in Fock space 



n(t)\n,t) = n\n,i). 



(50) 



Thus, in this space we can then write the partition func- 
tion as 



Z = Tr exp 



1 



2sinh(£fi) 



(51) 



In Ref. 0, [1[ one can define the occupation number for 
a frequency u>, Eq. (|42|) . Then by fitting the number of 
particles created as the usual Planck distribution Eq . (|4"3")) . 



one can then in principle fit the temperature of the radi 
ation. Here, we then choose to define (3 as 

91n(l + 1/N) 



P = 



duj 



(52) 



This implies that all of the time-dependence of the system 
is encoded into the temperature of the system. 

We can see that Eq. lfSTj) is just the standard entropy 
for a time-independent harmonic oscillator; however, the 
temperature here is time-dependent. Therefore the en- 
tropy is also time-dependent. 



B. Entropy 

Using Eq. (|5ip and the thermodynamic definition of en- 
tropy, we can then write the entropy of the system as 



S 



lnfl 



-/3u> \ 



-/3uj 



(53) 



Therefore, this is again just the entropy of the usual time- 
independent harmonic oscillator. From Eq.® it follows 
that the temperature is time-dependent. 

First we consider the entropy of the entire system. In 
Figure [T] we plot the entropy of the entire system as a 
function of dimensionless time t/R s . Figure [T] shows that 
the system starts with an initial entropy of zero. This is 
expected since initially there is only one degree of free- 
dom, meaning that S = m(l) = 0. Here we have normal- 
ized the initial entropy of the shell to be zero. To justify 
this normalization, consider a solar mass black hole. Un- 
der the usual Bekenstein-Hawking entropy, the order of 
magnitude estimate of the entropy of a solar mass black 
hole is Sbh ~ 10 75 . Now consider that the shell is actu- 
ally made up of protons. The initial entropy of the shell 
then is approximately Ss.o ~ 10 57 . Comparing the en- 
tropy of the final black hole versus the initial entropy of 
the shell, the entropy of the final black hole is much much 
greater than that of the initial entropy of the shell, thus 
the initial entropy of the shell only contributes a negligi- 
ble amount of entropy to the entropy of the final black 
hole. Thus our normalization of the initial entropy of the 
shell to zero is justified. As t/R s increases, initially the 
entropy increases rapidly, then settles down to increase 
approximately linearly. Due to the linear increase, we 
see that as t/R s goes to infinity, the entropy will then 
diverge. This is again expected since as the asymptotic 
time goes to infinity, the number of particles that are 
produced diverges (see Ref.Q)- This is a consequence 
of the fact that we keep the background fixed (i.e. R s 
is a constant). In reality, R s should decrease over time 
since the radiation is taking away mass and energy from 
the system. Therefore as t/R s goes to infinity, the en- 
tropy of the entire system as measured by the asymptotic 
observer diverges as R — ► R s . 

This is consistent with the results found in Ref. [13] . 
Here the authors consider the time-dependent non- 
equilibrium evolution of a black hole as well as the in- 
corporation of the given-off radiation. Here one can see 
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FIG. 1: We plot the entropy of the entire system as a function 
oit/R a . Except for the initial increase, the entropy increases 
approximately linearly as t/R s increases. 



that the entropy of the system diverges as the time goes 
to infinity. 

Now we consider just the radiation which is induced 
during the collapse. In Figure [5] we then plot the en- 
tropy as a function of t/R s . Figure [2] shows initially the 
entropy of the system is zero. Again, this is expected 
since initially the domain wall is in vacuum, meaning 
that initially there is no radiation being induced. There- 
fore the only degree of freedom is that of the domain wall, 
which then gives that the initial entropy of the radiation 
must be zero. As t/R s increases, initially there is a rapid 
increase in the entropy, but again, the entropy then in- 
creases linearly as t/R s increases. As in the case of the 
entire system, as t/R s goes to infinity, the entropy of 
the radiation also diverges. This is expected since the in- 
duced radiation diverges as R — > R s ; hence as the domain 
wall approaches the horizon, the occupation number for 
the induced radiation diverges. 
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FIG. 2: We plot the entropy of the induced radiation created 
during the collapse as a function of t/R s . Except for the 
initial increase, the entropy increases approximately linearly 
as t/R s increases. 



In Figure[3]we plot the entropy as a function of asymp- 
totic observer time t of both the entire system and the 
induced radiation during the time of collapse. Figure 
shows that except the initial increase in the entropy, for 
later t/R s (approximately t/R s = 7.5), the slopes of the 
entropy versus time are approximately equal. Therefore, 
one can expect that the entropy of the domain wall is 
approximately constant for late times. 
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FIG. 3: We plot the entropy as a function of t/R s for both 
the entire system and the induced radiation during the time 
of collapse. At the time t/R a « 7.5, the slope of the two lines 
is approximately equal. 

What is of interest is the entropy of the collapsing do- 
main wall, since this will form a black hole. To find the 
entropy of the domain wall, we can take the entropy of the 
entire system and subtract off the entropy of the induced 
radiation (since these are the only relevant objects which 
contribute to the entropy). The result is then given in 
Figure [4] We see some interesting features of the time- 
dependent entropy in Figure [H Here we discuss these 
features; first we will make some general comments on 
the overall behavior of the time-dependent entropy, fol- 
lowed by discussions of late-time (values of t/R s > 8) 
and the early-time (values of t/R s < 8) behaviors of the 
time-dependent entropy, respectively. 

Figure 2] shows that initially the entropy of the do- 
main wall is zero. As stated above, this is expected since 
initially there is only one degree of freedom. As t/R s in- 
creases, the entropy of the domain wall rapidly increases. 
However, for late times (t/R s > 8), the entropy of the 
domain wall goes to a constant. From the discussion 
above, this is expected since the late-time entropies for 
the entire system and for the radiation arc approximately 
parallel. This feature is somewhat artificial though: the 
entropy is constant since we arc assuming that the mass 
is approximately the Hamiltonian of the system, which 
is a constant of motion. This means that since we are 
holding the mass of the domain wall constant, we need 
to keep adding energy to the system to counteract the 
loss of mass (energy) from the radiation. Therefore one 
can expect that the entropy of the domain wall must be 
a constant for late times. 
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In reality, radiation takes mass away from the system, 
so the entropy of the domain wall will go to zero as R s 
goes to zero. Hence if the domain wall doesn't start off 
with enough mass, it could evaporate before ever creating 
the black hole. This means that after the domain wall 
disappears, all the entropy will go into the entropy of the 
radiation, since all the mass has dissipated. 



Entropy of Shell 




FIG. 4: We plot the entropy of the domain wall as a function 
of t/R 3 . Here we see for t/R s > 8 the entropy is constant in 
time with a value of S ~ 0.7 R 2 S . 

For large values of t/R s (t/R s > 8) the entropy of the 
domain wall is constant. From Figure 0] we see that our 
numerical value for the entropy of the domain wall for 
this region is 

S ps 0.7i^. 

Comparing with Eq.([T]), we can view this discrepancy as 
a shift in the Schwarzschild radius R s . In order to get 
the theoretical value for the entropy, Eq. (JTJ) , we see that 
we would require R s — > 2.11R S . This is an understand- 
able numerical error, which implies that our numerical 
solution is of the same order as the Hawking-Bekenstein 
entropy. 

For small values of t/R s {t/Rs < 8), the entropy of 
the domain wall seems to stop increasing around t/R s ps 
7.5. This means that the entire change in entropy of the 
domain wall occurs in the region < t/R s ps 7.5. At 
first this value for t/R s seems arbitrary, however there is 
some physical significance to this time. 

First, we notice that the volume of the domain wall 
becomes approximately constant at t/R s ps 7.5. To illus- 
trate this we will first require i?o = nR s , where n G R 
(R being the real numbers). From Eq. (fl6|) we can then 
write the position of the domain wall as 

R(t)=R s + ( n _ i) e -*/«.) . 

For illustration purposes we choose the value n = 10, 
which at time t/R s = 7.5 corresponds to a distance of 
R ps 1.005ii s with velocity R = -B ps -0.005 (the 



numerical values for the position and velocity will be 
even less for smaller values of n). Hence the value of 
t/R 8 ps 7.5 corresponds to the region where the domain 
wall is very close to the horizon and the velocity of the 
domain wall, with respect to the asymptotic, goes ap- 
proximately to zero. Therefore, as far as the asymptotic 
observer is concerned, the dynamical process of the do- 
main wall collapsing is over and the domain wall has 
come to rest, meaning that there is now approximately 
a constant volume. From this point on, classically it 
takes an infinite amount of time for the domain wall to 
reach the Schwarzschild radius, hence travel the distance 
R ps 0.005i? s , see Section HU 




FIG. 5: We plot 8 of the entire system and the induced radi- 
ation as a function of t/R s . Here we see for t/R a ps 7.5 the 
two values of /3 become approximately equal, meaning the two 
systems are in thermal equilibrium. 

Second, we can show that the entire system and the 
induced radiation come into thermal equilibrium at this 
time. In Figure [5] we plot (3 versus t/R s for the en- 
tire system (continuous curve) and the induced radia- 
tion (dashed curve) . Figure [5] shows that for the time 
t/R s ps 7.5 the values of the two /3's become approxi- 
mately equal, meaning that the entire system and the in- 
duced radiation are now at the same temperature. There- 
fore the system is now in thermal equilibrium, meaning 
that there is no more change in entropy of the domain 
wall as t/R s increases. Further more, the fluctuations 
(departure from thermality) in 8 become very small at 
this time, as discussed in Refs.@, H|. 

Finally we can evaluate the the chemical potential for 
both the entire system and for the induced radiation. 
From definition we can write the chemical potential as 



In Figure [6] we plot the chemical potential for both the 
entire system and for the induced radiation. We can 
see that as t/R s increases the chemical potential of the 
entire system and the induced radiation goes to zero. 
This means that the dispersion of particles goes to zero 
and the system goes into equilibrium. 
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FIG. 6: We plot fj, versus t/R 3 . The solid line corresponds 
to the entire system while the dashed line corresponds to the 
induced radiation only. Here we see that as t/R 3 increases, 
the chemical potential for each goes to zero. 

During the dynamical process, the entropy increases 
almost linearly. If one applies a best-fit line, we see that 
the entropy oscillates about the best line. These oscil- 
lations may be attributed to several different circum- 
stances. First, the oscillations may be caused by the 
non-thermal property of the radiation (see Ref.Q). Sec- 
ondly, these oscillations may be a manifestation of the 
error associated with the numerical calculations. Lastly, 
the oscillations may be an artifact of expanding the cal- 
culations beyond the region of validity, since we are using 
the near horizon approximation. Hence for values large 
compared to i? s , we cannot completely trust our result. 

V. CONCLUSION 

We studied the time evolution of the entropy of a col- 
lapsing spherically symmetric domain wall (representing 
a shell of matter) using the Functional Schrodinger and 
Louisville-von Neumann formalisms from the point of 
view of an asymptotic observer. This was done by inves- 
tigating the change in entropy for the entire system (the 
domain wall and the induced radiation) and that of only 
the radiation induced during the collapse semi-classically, 
i.e. the background of the collapsing domain wall was 
treated classically while the radiation was treated quan- 
tum mechanically. 

To summarize, we have arrived at the following re- 
sults. From Figured] we see that dS/dt > for the entire 
system, which is in agreement with the second law of 
Thermodynamics. Second, from Figure [2J we see that 
initially the entropy of the radiation increases at a lower 
rate than that of the entropy of the entire system; how- 
ever, for larger values of t/R s the increases of each are 
approximately the same (as seen in Figure [3]). Finally, 
from Figure |4] we demonstrated that the entropy of the 
domain wall is constant for late times. As discussed in 



Section IIV B[ this is expected since we are holding the 
mass of the domain wall constant, thus we must keep 
adding energy to the system to maintain the mass of the 
domain wall. In a completely realistic model, this would 
not be the case. The entropy would be expected to de- 
crease as R s decrease, since R s = 2GM. So as the mass 
decreases the Schwarzschild radius would also decrease, 
meaning that the entropy of the domain wall would go 
to zero as R s goes to zero. Here we demonstrated that 
the late-time entropy (large values of t/R s ) of the do- 
main wall does in fact go to a constant as predicted by 
Bekenstein, see Refs.[3l.[l3j. We found that the late time 
entropy of the domain wall is S « 0.7 'R 2 S . By comparing 
Eq.([T]) with our result, we see that our result is consistent 
with the accepted value for the entropy. 

Figure [4] displays some other interesting features. First 
we see that all the change in entropy occurs during small 
values of t/R s , with the increase in entropy ending at a 
time t/R s f=a 7.5. As discussed in Section HVB1 this time 
is not an arbitrary time. To see this, first, for different 
values i?o, this time corresponds to when the domain wall 
is very close to the horizon and the corresponding veloc- 
ity is very close to zero. From this point on the domain 
wall takes an infinite amount of time to traverse the re- 
maining distance to the horizon. Thus, as far as the 
asymptotic observer is concerned, the dynamics of the 
domain wall has finished by the time t/ R s « 7.5, hence 
the observer won't see any further change in the entropy 
since the volume of the shell is approximately constant. 
Second, the "inverse temperature" (f3) for both the entire 
system and the induced radiation become equal at this 
time. This means that the two systems come into ther- 
mal equilibrium, as seen in Figure Furthermore, the 
fluctuations in the temperature become irrelevant by this 
time, as discussed in 0, @|. Finally from Figure [6] we see 
that the chemical potential for both the entire system and 
the induced radiation go to zero near this time. Again, 
this signifies that the dispersion of the induced particles 
goes to zero and the system arrives at equilibrium. 

The increase during this region is approximately lin- 
early with oscillations about the best fit line. These os- 
cillations may be attributed to several different circum- 
stances. First, the non-thermal property of the radiation, 
as discussed in Ref.0]. Second, a manifestation of error 
associated with the numerical calculations. Lastly, we 
may have extended our calculation beyond our region of 
validity, since we are using the near-horizon approxima- 
tion, so we may not be able to completely trust the linear 
result in this region. 
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